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Abstract 

o 

Q 

Helicity amplitudes for the 'y*N ^ A transition are calculated using the 
I cloudy bag model. A correction for center-of-mass motion is carried out using 

T— I ! 

[ a modified Peierls-Thouless projection method. This reduces the magnitude 

[ of the transition amplitudes at small momentum transfer and enhances them 

I at modest momentum transfers. Our calculation shows that the pion cloud 



X 



contributes substantially to the transition helicity amplitudes, with the fi- 
nal result giving reasonable agreement with the corresponding experimental 
values. 

I. INTRODUCTION 

The nucleon-delta electromagnetic transition amplitude is an outstanding example of the 
success of the quark model |l|]. There have been many theoretical and experimental explo- 
rations of this transition process. In a naive quark model the 7*A^ — A transition occurs 
only by an Ml transition, while the E2 process is fully suppressed 0. In more sophisti- 
cated models, quarks can interact through, for example, one-gluon exchange in addition 
to the confinement potential between them. Then it is possible for configuration mixing, 
involving the excitation of one quark to a rf-state, to generate a small, but nonvanishing, E2 
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amplitude 0. To extract the 7*A^ A amplitude from experimental data is not an easy 
task. There are some uncertainties on the subtraction of background, and the results are 
somewhat model dependent |^. With the advent of new generation of accelerators, much 
more accurate measurements will be made. The prospect of the high quality data should 
test various hadron models well and help to build more realistic ones. 

The cloudy bag model (CBM) improves the MIT bag model |^ by introducing an 
elementary, perturbative pion field which couples to quarks in the bag in such a way that 
chiral symmetry is restored. The pion field significantly improves the predictions of the static 
properties of baryons. It also provides a convenient connection to the traditional language 
in nuclear physics. Kalbermann and Eisenberg (KE) first used the cloudy bag model 
to calculate the photoproduction of delta and Roper resonances at the threshold. They 
showed that the pionic corrections can substantially improve the predictions of the helicity 
amplitudes. Subsequently Bermuth et al. P] updated the KE calculation with some algebraic 
corrections and the appropriate normalization of baryon states. By comparing pseudoscalar 
(PS) and pseudovector (PV) coupling versions of the CBM, they found that the two are 
almost identical for the Ml transition, except for the seagull term, which actually spoils the 
good fit to the experimental data. The recoil effects were neglected in both calculations, 
which were limited to photoproduction at the peak of the delta resonance. 

In the standard treatment of the cloudy bag model, the cavity approximation for quark 
fields is adopted, and the pion field is then added in a perturbative way. The wave function 
of a bare baryon is just a direct product of individual quark wave functions, similar to 
the nuclear shell-model wave function (independent particle motion). This type of wave 
function is not a momentum eigenstate although the Hamiltonian commutes with the total 
momentum operator. The matrix elements evaluated between such static states contain 
spurious center of mass motion which ought to be removed. Early studies indicated that 
the correction for spurious center of mass motion is significant 0. It is expected to be 
most important in calculations where relatively large momentum transfers are involved. 
There are several intuitively motivated prescriptions P| for the correction of center of mass 
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motion, however, none of them are fully satisfactory. In this work, we have chosen to use 
the Peierls-Thouless (PT) method to eliminate the center of mass motion, since it is the 
most convenient for our purposes. It is basically a nonrelativistic approach and requires two 
consecutive projections of the product wave functions to form a Galilean invariant state for 
a composite baryon. The energy expectation value evaluated using the PT projected wave 
function has the correct nonrelativistic form, whereas this is not the case for Peierls-Yoccoz 
projection. 

In this paper, we calculate the nucleon-delta electromagnetic transition amplitudes with 
respect to the virtual photon. The spurious center of mass motion is corrected by using 
the PT projection method. As a first step, we assume exact SU (6) symmetry for the quark 
structure of the baryons, so that all quarks in the ground state of and A are in the 
s state. The notation of references |lTI| , p!2| is followed. We briefly review the method to 
construct the PT wave function in Sec. II. The calculation of helicity amplitudes for virtual 
photopro duct ion of the delta is performed in Sec. Ill, and in Sec. IV we present the numerical 
results. Finally in Sec. V we summarize our results. 



II. GALILEAN INVARIANT BARYON STATES 



We start with the chirally invariant Lagrangian density of the cloudy bag model 



C = {tq^df^q - B)ev - -qq^s 



(1) 



where Oy is a step function which is one inside the bag volume V and vanishes outside, and 
A5 is a surface delta function. In a perturbative treatment of the pion field, the quark wave 
function is not affected by the pion field and is given by the MIT bag solution 



g(r) 



ia ■ rf{r) 



e{R-r) 



(2) 
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where R is the spherical bag radius. For the ground state of a massless quark g{r) = 
Nsjo{uJsr/R), f{r) = NM^,r/R), where cu, = 2.0428 and A^^ = cu,/87ri?3j2(^^)(^^ _ 

The bare baryon is taken to be composed of three quarks with the spin-flavor wave 
function given by SU{Q) symmetry. The space component is naively the direct product of 
three quark wave functions in coordinate space 

^'(xi, X2, X3; x) = g(xi - x)g(x2 - x)g(x3 - x). (3) 

Here x indicates the location of the bag center, while Xi, X2, and X3 specify the positions 
of the three quarks. Clearly this wavefunction does not have definite momentum and is 
not a momentum eigenstate. A momentum eigenstate of the baryon can be constructed by 
making a linear superposition of the localized states, namely, 

^'py(xi, X2, X3; p) = N'{p) j d^xe*P-''^'(xi, X2, X3; x), (4) 

where the subscript PY stands for Peierls-Yoccoz projection |]13|, and N'{p) is a momen- 
tum dependent normalization constant. It can be shown that \1/py(p) = e*P'''=m\l/in(p), where 
Xcm = (xi + X2 + X3)/3 is the the center of mass of the baryon and ^&m(p) is the appropriately 
defined intrinsic part of the wave function. Since \I'm(p) still depends on the cm. momen- 
tum, it violates translational invariance. To overcome this problem, Peierls and Thouless 
(PT) 0] proposed to make another superposition of these momentum eigenstates, i.e., 

vI/pT(xi,X2,X3;p) = Nip) y'c/V«^(p')e^(P-P')-''^-vl>pY(xi,X2,X3;p'). (5) 

The weight function, w{p'), should in fact be chosen to minimize the total energy, but this 
would be quite complicated to implement here. Instead, we choose w{p') = 1 for simplicity 
and convenience. Then integrations over x and p' can be carried out easily. This leads to a 
much simplified PT wave function, 

^fpTlxi, X2, X3; p) = A^pTe'P"''^-g(Xi - Xcm)g(x2 - X„n)g(x3 - Xcm), (6) 

where the normalization factor, iVpT, is given by the condition 
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1/2 

(8) 



yc/^xic/^a;2c/^a;3^J.T(xi,X2,X3;p')*PT(xi,X2,X3;p) = {27rf6^^\p-p). (7) 
This leads to 

= [3 l^r,^rMrM-r. - r, 

where p(r), the quark density at the position r, is defined as 

p(r) = q\r)q{v) = [g\r) + f{r)] e{R - r). (9) 

Notice that, for this simple version of the PT projection, Npj^ is a momentum independent 
constant and the wavefunction in Eq. (|^) is manifestly Galilean invariant. 

III. THE HELICITY AMPLITUDES IN THE CLOUDY BAG MODEL 

From the CBM Lagrangian given in Eq. (||), the conserved local electromagnetic current 
can be derived using the principle of minimal coupling dfj_ + iqA^, where q is the charge 
carried by the field upon which the derivative operator acts. The total electromagnetic 
current is then 

J^(a;)=j,^(x)+j^(x), (10) 
j!:ix) = J2Q,eUx)rqaix), (11) 

a 

j^{x) = -ie[TT\x)d''TT{x) - 7r(x)9M(x)], (12) 

where qa{x) is the quark field operator for fiavor a, Qa is its charge in units of e, and e = |e| 
is the magnitude of the electron charge. The charged pion field operator is defined as 

7r(a;) = -^[7ri(x) + i7r2{x)], (13) 

where n{x) either destroys a negatively charged pion or creates a positively charged one. 

It is customary to define the helicity amplitudes for the electroproduction of the delta 
as [HI 
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^3/2 = -^{A;sA = 3/2| J(0) ■ e\N-SN = 1/2), (14) 

A/2 = ^^(A; SA = 1/2| J(0) ■ e1iV; sn = -1/2), (15) 

where the A is at rest and the photon is travelhng along the z-axis with right-handed 
polarization, e = h 0). The spin projections of A and along the z-axis are denoted 

as sa and sn respectively. For a virtual photon, the three- momentum in the A rest frame 
is given by 



(Mi -M^- Q2^2 
4M? 



with Q = \l~—q^ the magnitude of the four momentum transfer. The photon energy is related 
to this by go = l'?!^ ~ Q^; where for a real photon we have = 0, so that uj^ = \qo\ = 
{Ml - MI)/2Ma. 

The experimentally extracted, resonant, helicity amplitudes are to be associated with the 
fully dressed initial and final baryons. In the cloudy bag model, due to the ttBB' coupling, 
a physical baryon state is described as a mixture of a bare bag and its surrounding pion 
cloud, 

\A) = ^/zi[l + (Ea ~Ho- AHjA)-'Hj] \Ao), (17) 

where is the bare baryon probability in the physical baryon states, A is a projection 
operator which projects out all the components of \A) with at least one pion, and Hj is the 
interaction Hamiltonian which describes the process of emission and absorption of pions. 
The matrix elements of Hj between the bare baryon states and their properties are then 



given by |[TT,12 



v^fik) ^ {Ao\H,\n,{k)Bo) = ^ p^;^^^EC^S^(S:;, ■ k)CSl^-J: ■ e,), (18) 
w^fik) ^ {Ao7r,{k)\Hj\Bo) = [<^(^)]* = -v^f{k) = v^f{-k), (19) 
where the pion has momentum k and isospin projection j, f^^ is the reduced matrix element 



for the ttBo — > Aq transition vertex, u{kR) = 3ji{kR)/kR, ook = \J k'^ + and Sm and U, 
are spherical unit vectors for spin and isospin, respectively. 



Under the approximation that there is at most one pion in the air, there are three different 
processes contributing to the 7*A^ ^ A vertex, as shown in Fig. 1. Substituting the quark 
current operator, Eq. (11), and the physical baryon states, Eq. (17), into Eqs. (|T^) and (f5|), 
we have 



(A, saIj, ■ e|A^, sn) = yZ^Z^ [(Aq, sa| jg ■ e |A^o, sat) 

+ (Ao, sa\Hj{E -Ho- AHjA)-% ■ e{E -Ho- AHjA^^Hj \No, sn) 



(20) 



where E is the total energy of this transition process. The first term, as illustrated in Fig. 
1(a), is the quark core contribution and the second term, corresponding to Fig. 1(b), is the 
contribution from the 'yqq interaction with one vr in the air. By inserting two complete sets 
of states before and after the operator jq in the second term and using the PT wavefunctions 
given in Eq. (H), we obtain 



(/ 



NN\2 



dkk^u^kR) 



2/25 



5/4 



+ 



+ 



1 



(21) 



(22) 



where 5 = m/\ — itln, and is the renormalized ttNN coupling constant. The four terms 
in the right-hand side of Eq. (|^) correspond to four possible intermediate states, (iVA), 
(AA), (AA^), and (NN), respectively. The recoil corrected bare 'jNo —* Aq transition 
amplitude is 

e J^drr^g{r)f{r)ji{qr)K{r) 



^bare i^Q 



(23) 



TTi/BZoj; J^drr'^p{r)K{r) 
where K{r) = J d^xp{x)p{—x — f) is the recoil function to account for the correlation 
of the two spectator quarks. The renormalization constants, Z^, are determined by the 
normalization condition for the physical baryon state, i.e. 

d 



(24) 



E=mA 



where S"^, the self energy of the baryon A, is given by 



where P indicates that the principle part of the integral is to be taken. In this work, we have 
adopted the usual philosophy for the renormalization in the CBM. Throughout this work 



approximate relation, ~ \jt^) always used. There are uncertain corrections 

on the bare coupling contant f^^ i such as the nonzero quark mass and correction of center 
of mass motion. Therefore, we use the renormalized coupling constant in our calculation, 
jNN ^ 3 which correspond to the usual ttA^'A'" couphng constant, /^tvat — 0.081. As 



a result, the factor ^ Z2 is absorbed into the renormalized coupling constants in Fig. 
1(b). This treatment is equivalent to the original CBM formalism up to order {f^^)^ and 
consistent with current conservation. 

To evaluate the contribution caused by the photon-pion-pion coupling vertex [see Fig. 
1(c)], we use the usual plane wave expansion for the quantized pion field 

nj{x, t^O)^J ^(^^^f^^^p/, [a,{k)e^>^-^ + a]{k)e-^'^-^] , (26) 

where aj{k) {aj{k)) annihilates (creates) a pion with momentum k and isospin j. The pion 
current operator becomes 



33 



(27) 



The transition amplitude at position x is just the matrix element oi Jt^^x) evaluated between 
the physical baryon states. Using the identities 

a^{k) \A) = {Ea -Uk- H)-%{k) \A), (28) 
aAk')a,{k) \A) = {Ea - uj^ - - H)-'V^,{k){EA - Uk' - H)-%,{k') 

+ {Ea -cok- cok' - H)-%,{k'){EA - co^ - H)-%{k) \A), (29) 

we have 
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X 



Y: [vf'jik', k)G''{k', k) + ryj^.,(fc, k')G^{k, k') 



(30) 



Here B denotes the intermediate baryon states (restricted to and A here), and 
f^BfNB u{kR)u{k'R) 



VrAk',k) 



G^ik', k) 



ml 167r3(a;fca;fc/)V2 

tB 

1 1 1 



■3)1 



SB 



G^{k\k) = 



1 1 



+ 



{uJk + uJk' + 6){uJk + S) {oJk> + 6 - uj^) {ujk + 
1 

+ 



(31) 
(32) 

(33) 



{Uk' + 6 — UJ^){uJk + OJk' — OJ^) 

G'^{k, k') and G^{k, k') are obtained by the interchange of k and k' in the corresponding 
equation. The three terms appearing in Eqs. (|32D and ( ^3]) correspond to the three different 
time orders in the time-ordered perturbation theory, as illustrated in Fig. 1(c). Using the 
translational invariance of the electromagnetic current operator, j^{x) = e*^'^j'^(0)e~*^'^, 
then the 7*iV — > A helicity amplitudes due to the •jim interaction are simply given by 



A(Q2) = I £xe''^-'^{A, sa\L{^) ■ e\N, sn)- 



(34) 



After performing some spin and isospin algebra, we obtain 

(/^^)'|g1 fd^kk^sm^eu{kR)u{k'R) 



^3/2 (Q^) 



240V6^vr3m2 
G^{k, k') + 3G^{k', k) + 2G^(fc, k') 

U^^f\^ fd^kk^sm^eu{kR)u{k'R) 



UJkUJk' 



UJkUJk' 



X 



2G^{k', k) - G^'ik, k') + G^(A;', k) + 4G^(A;, k') 



(35) 



(36) 



where k' = k + q, uJk' = y k'"^ + m^, and 6 denotes the angle between k and q. It is worthwhile 
to mention that the form of our results for Fig. 1(c) are expressed in a quite different form 
than those of KE and Bermuth et al. p[ where the integral variables are k and k' in 
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their formulation. We believe that our expressions are more straightforward and manifestly 
respect the three momentum conservation at the '-firn vertex. Our numerical results for this 
contribution also appear to differ from those of Refs. . 

IV. NUMERICAL RESULTS 

Figure 2 shows the recoil function K{r) due to the two quark spectators. It cuts down 
the contributions from the quark wave functions near the bag boundary by up to 50 %. The 
overall effect of this recoil correction on the bare bag contribution to the typical '~f*N —>■ A 
helicity amplitude, A3/2, is shown in Fig. ^. In the real photon limit {Q^ 0), the magnitude 
of the 7*A^ —>■ A transition amplitude increases with the bag radii in a fashion similar to 
that of the magnetic moment of bare baryons. The correction of the center of mass motion 
usually reduces the bare transition amplitudes by 5 to 10 % for ^ 0.5 GeV^ within a 
reasonable range of bag radii. However, this recoil correction would flip sign and increase 
the transition amplitude for larger momentum transfers. 

The expressions of A^''^ and A^'^^ (from Fig. 1(b) and 1(c) respectively) involve pion-loop 
integrals, which contain poles originating from the baryon mass difference in the propagators. 
Using the relation P /g"^ = 0, it can be shown that 

/-"^ dkfjk) f-- fik)-f{ko) tiT 

}^Jo k^-kl-re =L k^-kl ^2kJ^^'^^ 

where fco > 0. To be consistent with the phenomenology of the cloudy bag model, we limit 
our calculation to the small momentum transfer region (approximately ^ 1.6 GeV^). 
For higher momentum transfers, relativistic effects are essential and the support of the 9 
integral is considerably more complicated. Consequently, for our purposes, there is only one 
relevant pole, ojk = cu^, for Fig. 1(b) and an extra pole, ujk' = 1^7, for Fig. 1(c). 

Fig. ^ shows the individual contributions to the helicity amplitude, A3/2, from the various 
diagrams for a typical bag radius, R = 0.8 fm. It is clear that the contributions from the 
pion cloud are signiflcant in the cloudy bag model. As — > 0, the pionic effects account 
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for two thirds of the total amphtude, or roughly twice that of the quark core for the chosen 
bag radius. As increases beyond 0.2 GeV^, the contribution from the 77r7r interaction 
decreases rather rapidly, so that it has nearly vanished at ~ 1.2 GeV^. For Q"^ ^ \ GeV^, 
the pionic effects contribute about 40 % of the total amplitude. Notice that, since the A is 
a resonance, the helicity amplitudes are actually complex. In the CBM the imaginary part 
arises from thresholds in the pion loop integral, and amounts to 15 ~ 20% of the real part. 

The real parts of total helicity amplitudes, A = A^*^-* + A*-^-* + A^*^-*, are presented in Fig. |^. 
With the contributions of the pion cloud, the bag radius dependence is quite different from 
that for the bare transition amplitudes shown in Fig. 3. This can be explained by the 
fact that the pionic contribution is competing with that of quark core, since a small bag 
radius means a strong pion cloud. In the small region, the pion cloud compensates more 
than the loss in the bare amplitude when using a small bag radius. Generally, the smaller 
the bag radius, the larger the total transition amplitude. We list the helicity amplitudes 
corresponding to the real photon limit at the A resonance in Table 1. With the small bag 
radius i? = 0.7 fm, we are able to reproduce the experimental helicity amplitude in this 
model. 

Finally, in Fig. |^, we show the real part of the invariant 7*A^ ^ A magnetic transition 
form factor, G^{Q'^), in comparison with the experimental measurements [T^. The relation 



between the helicity amplitudes and the invariant transition form factors are given in Ref. 



| r?| . With the bag radius R = 1 fm, our calculation agrees reasonably well with the data 



for modest Q^, but is less satisfactory at — GeV^ and in the region ^ 1 GeV^. At 



the real photon limit, we are able to get the PDG [T^ value with R = 0.7 fm, however, at 
nonzero the prediction with this bag radius would not be in agreement with the data. 
The center of mass correction reduces G^^(O) by approximately 5 % and makes the Ml form 
factor somewhat harder, similar to what was found in the case of the nucleon magnetic form 
factor [|1|]. 

In the present work the quark cores of nucleon and delta were assumed to have SU(6) 
symmetry, i.e., there is no deformation assumed, so that the spatial parts of the wave 
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functions for the and A are purely s state. Hence, here the only possible source for an 
E2 amplitude is from Fig. 1(c). Using Eqs. (21), (22), (35) and (36), we obtain E2/M1 = 

(Ai/2 - A3/2/V3)/(Ai/2 + V3A3/2) 0.3 X 10"3. Because of the severe cancellation in the 

numerator, this result is about two orders of magnitude smaller than the phenomenological 
approach P|,p!5[ (The latest estimate by Particle Data Group is —0.015 ± 0.004). 



V. SUMMARY 

In conclusion, we have calculated the 7*A^ A transition form factors in the cloudy bag 
model, including the center of mass correction via a simplified Peierls-Thouless projection 
method. The effect of this recoil correction is to slightly reduce the magnetic form factor at 
small momentum transfer and to make the form factor slightly harder. Generally, with the 
FT projection the transition moment is reduced by about 5 ~ 8%. 

The pion cloud contribution proved to be crucial to account for the measured helicity 
amplitudes using a reasonable bag radius in this model. In similar calculations using con- 
stituent quark models (the nonrelativistic [0] and relativized quark models |1^), the helicity 
amplitudes are usually significantly underpredicted. 

Because of the large range of momentum transfer measured, we have been unable to 
obtain a totally satisfactory description of all of the data using our essentially nonrelativistic 
formalism. Typically the form factors are too stiff. Unfortunately, there is as yet no covariant 
formalism for quark bag models or soliton bag models. Using a relatively simple prescription 
for Lorentz boosts [jl9|, the invariant form factor G^{Q'^) has shown improved behaviour. 



Obviously a crucial next step will be to construct a quark model with improved Lorentz 
transformation properties, in addition to the inclusion of a perturbative pion cloud. It will 
also be important to explore the effect of SU{Q) violating admixtures in the baryon wave 
functions 0. 

This work was supported by the Australian Research Council. 
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TABLES 

TABLE I. Helicity amplitude of delta photoproduction, A3/2, in units of 10~^GeV~^/^. Here 
static denotes the static calculation and PT denotes the Peierls-Thouless projection. The indices 
a, b, and c correspond to the Figs. 1(a), 1(b), and 1(c) respectively. The latest estimate by Particle 
Data Group is —258 it 6. 





static 


PT 


R(fm) 


a 


b 


c 


total 


a 


b 


c 


total 


1.0 


-115 


-53 -21i 


-64 -18i 


-205 


-106 


-49 -19i 


-64 -18i 


-195 


0.9 


-106 


-62 -20i 


-80 -20i 


-216 


-98 


-58 -18i 


-80 -20i 


-205 


0.8 


-96 


-73 -19i 


-101 -21i 


-233 


-86 


-68 -17i 


-101 -21i 


-222 


0.7 


-85 


-87 -17i 


-129 -22i 


-260 


-79 


-80 -16i 


-129 -22i 


-249 
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FIGURES 

FIG. 1. Diagrams illustrating the various contributions included in the calculation. The inter- 
mediate baryons B and B' are restricted to the N and A here. 

FIG. 2. The recoil function K(r), defined following Eq. (23), due to the two quark spectators 
in arbitrary units. 

FIG. 3. The efi'ect of the center of mass correction on the helicity amplitude, ^3/2) for the bare 
bag. The number on each curve indicates the bag radius in fm for the calculation. 

FIG. 4. Individual components for the helicity amplitude, ^3/2, with a typical bag radius R 
= 0.8 fm. 

FIG. 5. The real parts of total -f*N A hehcity amplitudes, Re[Ay2{Q'^)] and Re[Ai/2{Q'^)]. 

FIG. 6. The real part of 7* A/" A invariant Ml transition form factor, Re[G^{Q'^)]. Four 
different bag radii are used in the calculation. The experimental data are taken from Refs. [16,17]. 
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